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We extend the notion of Dirac oscillator in two dimensions, to construct a set of potentials. 
These potentials becomes exactly and quasi-exactly solvable potentials of non-relativistic quantum 
mechanics when they are transformed into a Schrodinger-like equation. For the exactly solvable 
potentials, eigenvalues are calculated and eigenfunctions are given by confluent hypergeometric 
\ functions. It is shown that, our formulation also leads to the study of those potentials in the 

. framework of the supersymmetric quantum mechanics. 
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t^J- 1 INTRODUCTION 

(N : 

The solution of the (2 + l)-dimensional Dirac equation are of special interest, because of the rapid growth in 
nanofabrication technology that has made possible to confine laterally two-dimensional (2D) electron systems. These 
quantum confined electron systems are referred to as artificial atoms where the potential of the nucleus, in the non 
relativistic case, is replaced by an effective potential of the form V — \r 2 which is often used as realistic approximation 
0,0 El- Although the parabolic potential appears to be a good approximation for artificial atom structures, their 
modelling with various potential profiles will be interesting from the theoretical point of view as well as from its 
, practical applications. 

In order to analyze relativistic effects on the spectrum of such physical systems one should construct Dirac equation 
including adequate potentials and obtain its solution. For the relativistic case, the spectrum and properties of such 
systems can be determined by using two dimensional Dirac oscillator^ 0, Relativistic extensions of the various 
exactly and quasi-ex actly solvable (QES) potentials have also turned out to be of importance in the description of 2D 
phenomena |a YA lit M. 1 1 0| . Different condensed matter physics phenomena point to the existence of (2 + l)-dimensional 
systems whose spectrum determined by Dirac equation Hamiltonian including various potentials. It is well known 
that the dirac equation is used for the description of spin-1/2 relativistic particle. Meanwhile we mention here that 
the Hamiltonian in the form of the Klein-Gordon equation, so called Feshbach-Villars equation, has been constructed 
in a two component form for spinless particles and in an eight-component form for spin-1/2 particles [l2l US . fl4 | . 
Regrettably, the Dirac equation is exactly solvable only in a very restricted potentials. It is the purpose of the present 
article to construct a Dirac equation including a class of potentials whose spectrum can exactly be determined. For 
this purpose we transform the Dirac equation into two Schrodinger-like equation, because there exist a large number 
of papers discussing the 2D, few electron systems, most of which are tackling the problem in the framework of the 
Schrodinger-like equation. 

A Dirac equation with an interaction linear in coordinates was considered long ago|l5j and recently rediscovered 
in the context of the relativistic many body theories 16]. The equation is named Dirac oscillator, since in the 
non-relativistic limit it becomes a harmonic oscillator with a very strong spin-orbit coupling term. Dirac oscillator 
has attracted much attention and the concept gave rise to a large number of papers concerned with its different 
aspects [171 Hil flit I20I I21I I22I I23L l24j . Analogous to the Dirac equation, with a modified momentum operator, which 
in the nonrclativistic limit turns out be the usual Schrodinger equation. As we have already noted the Dirac equation 
including various potentials might attract much attention because it may have some physical applications, particularly 
in the condensed matter physics. It seems that one can present more realistic models for the artificial atoms using 
the procedure given here. 

In order to solve Dirac equation, in this paper, we use functional approach which have been applied to solve 
Schrodinger equation for a exactly or QES potential profile. For a QES potential it is possible to determine alge- 



braically a part of spectrum but not whole spectrum |25l 1261 127| . Our approach also gives a hint to the solution of 
the problem in the framework of the SUSYQM. 

The method presented here consist of the followings. In section 2, we introduce (2 + l)-dimcnsional Dirac equation. 
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Using the structures worked for the Dirac oscillator we develop a method to construct a class of exactly and QES 
potential profile. We transform the Dirac equation into the form of the Schrodinger equation. In section 3, we 
construct a Dirac equation including, harmonic oscillator, Coulomb and Morse potentials. We obtain the corresponding 
eigenvalues and eigenfunctions. In section 4 a class of QES potentials are constructed and their ground state wave 
functions are explicitly determined. We conclude our results in section 5. 



METHOD 



The (2 + 1)— dimensional Dirac equation for free particle of mass m in terms of two-component spinors ip, can be 
written as 



c/3jipi+/3r 



(1) 



Since we are using only two component spinors, the matrices (3 and (3ji are conveniently defined in terms of the Pauli 
spin matrices which satisfy the relation <z,<7j = <5,y + ejjfeO/., given by 



/37l=<7i; /?72 = cr 2 ; /3 = er 3 . 



(2) 



In (2 + 1)— dimensions, the momentum operator pi is two component differential operator, p = — ih(d x , d y ), for free 
particle. In the presence of the magnetic field it is replaced by p — * p— eA, where A is the vector potential, and 
the 2D Dirac oscillator can be constructed by changing the momentum p — > p— imua^rr. We are now seeking for 
a certain form of the momentum operator that can be interpreted as exactly solvable Schrodinger equation in the 
non-relativistic limit. For this purpose we introduce the following momentum operator 



p — > p— eA+i<T 3 v(r)r . 

Then Dirac equation JQ) with the momentum operator J2J takes the form 

[E - C7 mc 2 ] ip = ccr + [p x + ip y - e (A x + iA y ) —i(v x + iv y )} ifj + 
cer_ [p x - ip y - e (A x - iA y ) + i(v x - iv y )] ip 



(3) 



(4) 



where pi,Ai, and Vi are the i th components of the momentum, vector potential and the potential in Cartesian 
coordinate system, respectively. In polar coordinate, x = r cos <f>, y — r sin 0, with the choices of the vector potential 
A. x = — ^4(r)sin0, A y — A(r) cos (j> and the potential v x = v{r) cos (f>, v y = v(r) sin cj) the 2D Dirac equation (0J 
takes the form 



E-ip+ = ce 



ce 



-i<p 



+ _ A eA \n + v(r)) ip- 
or r 0(p 

■ t- & h 9 ., ., , , 

-ih- — + i{eA{r) + v(r)) fa 

or r 06 



(5a) 
(5b) 



where E± — E ± mc 2 and tp± = ip±{r, <p) are the upper and lower components of the spinor ip. The substitution of 
the wave functions 



V>±(r, 



-f±(r) 



leads to the following set of coupled differential equations 

E-f+(r) = 



f d hi \ 

-icih— + y + (eA(r) + v(r)) j /_ (r) 
E+f-(r) = -ic(h^-^-(eA(r)+v(r))^f + (r). 

Our task is now to transform (|7a|) and l|7b|l in the form of the Schrodinger-like equation. Substitution of 

I {eA{r)+v(r)) 



W{r) 



E± = ±iche±; e = £+£-, 



(6) 

(7a) 
(7b) 

(8) 
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into l|7a(l and JTb)), leads to the following expressions 



e_/+(r) - (|: + W(r))/-(r) (9a) 



= f-|; + W(r))/ + (r). 



(9b) 



Notice that the result l|9a|) and (|9b|) which seem to hint at a supersymmetric treatment of the Dirac equation, because 
the supersymmetric operators can be expressed as A^ 1 = (t^: + W(r)) . It is obvious that the functional form of 
the superpotential W(r) and v(r) are the same except that the radial function ^. Thus the superpotential of the 
non-rclativistic quantum mechanics can be recognized as the potential of the rclativistic quantum mechanics. It is 
obvious that the expressions i|9a|) and Ij9bf> can be written in the form of the Schrodinger-like equation, by eliminating 
f+(r) and/or f-(r) between (I9all and 19b|l provides the following expressions 



(-£-^ 2 W + ^'(r)+e 2 ^/_(r) = 



(10a) 



-^ + W 2 (r)+W\r)-e 2 ^U(r) = 0. (10b) 

This is indeed an interesting result for the potentials V± = W 2 (r) ± W'(r) leading to a common spectrum, thus 
forming an isospectral system. 

In the following we analyze the solutions and the energy spectrum of the (2+l)-dimensional Dirac equation including 
various potentials. 

EXACTLY SOLVABLE POTENTIALS 

Appropriate choices of the superpotential W(r) permits the construction of the equations having exactly solvable 
potentials. In this section we illustrate that for a large class of potentials, the Dirac equation in the form of Schrodinger 
type equations of I jlOajl and ijlObJl possess exact solutions. 

Harmonic oscillator 

Let us start with the well known problem, namely Dirac oscillator. The Dirac oscillator can be constructed with 
the choices of the superpotential W(r) = ^ojxr — sti. Thus, the Schrodinger type equations l|10a|) and (|10b(l takes 
the form: 

(-^ + ^ + (^ T )V + ^ T (« + 3 )-e»)/_M = (lla) 

(-^ + ^^ + (>)V-^ + 1) -/)/ +W - 0. (Ub) 

In this case the vector potential A(r) and scalar potential v(r) are given by 

A{r) = \Br, w (r) = ^wr- 2£±I, (12) 
2 n r 

and the frequency lot can be expressed in terms of the Larmor frequency as follows 

w T =w + w £ =w+- — . (13) 
2m 

In order to solve (|llap . we change the variable z — ^ustt 2 , and introduce the wave function 

/_(*) = Cz^e-* g-{z) (14) 
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where C is the normalization constant, then 111 all takes the form 



r d 2 










+ n 









9-(z) = 0. 



The natural number n and energy satisfy the relation 

72 



9 E — m c to 



(15) 



(16) 



Note that the non-relativistic limit of the energy is obtained by setting E = E nr + mc 2 and considering E nr <C rac 2 , 
we obtain the non relativistic energy 

E nr = AnfkVT- 

We now investigate the dependence of the energy on the spin. In order to analyze spin effects, it is worth to obtain 
non-relativistic form of the l|llaf> : 



tf_cP_ e{£+l)h 2 1 

2™ dr 2 2mr 2 2 T T 



Enr /-(r)=0. 



(17) 



The corresponding Hamiltonian is the Hamiltonian of a harmonic oscillator with an additional spin dependent term, 
hu T (t + |) ■ 

To complete our analysis we turn our attention to the normalization of the wave function. It is easy to see that the 
solution of l|15|) is the associated Laguerre polynomials, L n 2 (z), then / (z) can be written as 



f-(z) =Cz^e-^L n +l2 (z). 
The upper component, f+(z), of spinor tp(r), can be obtained from the relation (|9a|l and it is given by 



(18) 



f+(z) = ^— z 2 e ?L n _\(z). 

Normalization condition in polar coordinate is given by 

(^)=^°°(l/ + (r-)| 2 + |/-(r)| 2 )dr=l. 
Associated Laguerre polynomials satisfy the orthogonality condition: 

J o z e z L e n (z)Ll(z)dr = j 

Thus, we hnally obtain an expression for the normalization constant C: 

^/UJ^elT^ + 1) 



(19a) 



C 



(tilot + elT(n)) T (n + I + f ) 



(20) 



We have solved Dirac oscillator in two dimensional space which leads to a series of interesting results. The Dirac 
oscillator has various physical applications particularly in semiconductor physics [2flj . 



Coulomb Potential 



Another well known examples of the exactly solvable Dirac equation is that relativistic Hydrogen atom. The problem 
can be solved exactly when the magnetic field is zero, A(r) = and in the presence of the Coulomb interaction |28ll34| 
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v(r) = 



h ( 2l+1 } \ n this case the superpotential is W(r) — 4^. eo ™l^ +1 ) — ^r-, and then the Schrodinger-like 



equations l(l()a[l and (|10b(l takes the form 
/ d 2 £(£ + 1) 



dr 2 



+ 



me 



+ 



me 



2ne h 2 r \Aire h 2 {£ + 1) 



c?_ (£ + !)(£ + 2) e 2 
dr 2 r 2 r ' \4ire H 2 (£ + 1) 



-e 2 /_(r) = 



-£ 2 /+(r) = 



(21a) 
(21b) 



Following the similar developments used in the construction of the harmonic oscillator problem, the equation l|21a|) 
can be transformed in the form: 



d 2 

z— + (2£+2~z) + n 
az z 



g-(z) = 



by changing the variable r = 2ir£ h 2 (n + £ + \)z/me 2 and the wave function 

,U(z)=Cz e+1 e-i g.(z). 
Natural number n and energy of the Hamiltonian satisfy the relation 

o E 2 — m 2 c 4 



h 2 c 



Ane h 2 (£ + I) J \ine h 2 {n + £ + 1) 



Solution of the 1)22(1 leads to the following expression for the wave function 

.U(z) = Cz^e-Ulf+\z) 

and from the relation l|9a[) we obtain 

f + (z) = =j**V* [ 2 (/+ l)L^{z) + (£ + 1 

Using the identity 



27T£ ?l 2 



(2£+l))L^(z) 



z^e-*L^{z)L f n {z)dr = (2n + t + lf +2 T ( n + ? + l \ 



after some straight forward calculation we obtain the normalization constant 

"2vre ft 2 {4(£+ l) 2 



C = 



me 



V e 2 _T{r 



(2n + 2£ + 2y l+i K T{n + 2£ + 2 



where K is given by 



K= 1 



me 
2ne h 2 



(2£ + l] 



(22) 



(23) 



(24) 



(25a) 



(26) 



(27) 



(28) 



(29) 



Due to the recent interest in the 2D field theory in the condensed matter physics, the 2D Coulomb potential is 
physically relevant and the results obtained in 2D exhibit some new features [6J, |3Jj . 



Morse Potential 



The Morse oscillator is exactly solvable quantum mechanical problem and it is used to model the interaction of 
the atoms in the diatomic molecules. In order to obtain its relativistic form we apply the standard procedure. The 
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choices of the parameters v(r) 



hi 



hae 



hb, A(r) — and W(r) — b~ae ar , provides the following potential 



b 2 -e 2 /_(r) = 



--^L+a 2 e- 2ar -a(a + 2b)e- c 
or 1 

^ • <r, 2 "' ha (a 2b) t /.(,■) = n 



In order to solve (|30a[) we change the variable e Qr = and then introduce the wave function 

f_(r) = z^- n e-i g-(z) 



then we obtain 



*'dz 2 



2b 

2n + 1 - z 

a 



g-(z) = 0. 



Energy expression for the Morse oscillator is given by 



e 2 = an(an — 2b). 

The wave functions can be obtained from the equations l|32l) and (|9a|l and are given by 

-In, 



m = 

f+{z) = z 



2 a 



nLr, 



, , 2&_ 2n +l, s 

(z) + zL£_ x (z) 



Normalization condition yields the following expression for the normalization constant 

e 2 _T(n)T(n- 1) 



C 



r (2k + 1 - n ) (e 2 _T(n) + 2a 2 T(n + 1)) 



(30a) 
(30b) 

(31) 
(32) 



(33a) 
(33b) 



(34) 



The potentials we have derived here is significant from both physical and mathematical points of view. For instance 
the relativistic quark model requires the solution of the dirac equation containing single quark potential. The potential 
behaving like r at a large distances and 1/r at a short distances has been recently treated by Muci [3(| • We mention 
that the Morse oscillator potential can also be used as a large distance potential in some various physical systems. 
Consequently we have constructed relativistic version of the three well known potentials whose eigenfunctions are 
associated to the confluent hypergeometric functions. In the following section we construct QES potentials. 



QES POTENTIALS 

The formulation given in section 2 is also useful to construct the relativistic version of the QES potentials. In 
this section we construct three of them. At this stage we mention that the underlying idea behind the quasi-exact 
solvability is the existence of a hidden algebraic structure. Our task is now to demonstrate by appropriate choices of 
relativistic potential we can construct QES Dirac equations. Our examples includes anharmonic oscillator potential, 
radial sextic oscillator potential and perturbed Coulomb potential. 



Anharmonic Oscillator potential 

The anharmonic oscillator potential have been widely used in many physical and chemical applications. In order 
to construct anharmonic oscillator potential we introduce 

M 

v(r) = h hur + hbr 2 + ha; A(r) = KBr; W(r) =a + uj T r + cr 2 . (35) 

r 



With these choices we obtain the following Schr6dingcr-like equations 

d 2 



. )j 2 , V++cu T + a 2 +e 2 ) f+(r) = (3(ia) 
+ V_ -Lu T + a 2 +s 2 I /_(r) = (36b) 



<9r 2 

where V± are given by 

V± = 2{au T ± 6)r + (2a6 + uj 2 t )t 2 + 2buj T r 3 + b 2 r 4 (37) 
and its groundstate wave function is given by 

/i 0) (r) = cxp (- \br 3 - \uj T r 2 - ar) ; (38) 



3 2 

A part of the spectrum of the anharmonic oscillator potential can be obtained in the framework of the QES problem 
or its approximate solution can be obtained by using perturbation theory. 

Sextic Oscillator potential 

Another well known QES potential is the sextic oscillator potential. Its relativistic form can be obtained by the 
following choices of functions: 

2M £ 

v(r) = h tvjjr + hbr 3 ;A(r) = hBr; W(r) = h cu T r + br 3 . (39) 

r r 

In this case we obtain the following equations: 

° 2 + ^±^ + V + +Af + (r) = (40a) 



92 ' + V- + e 2 ) /_(r) = (40b) 



dr 2 
where 

V± = {uJt ± b(2£ T 3))r 2 + 2bu T r A + b 2 r 6 - io T (2£ =p 1) (41) 
Its groundstate wave function is given by 

/i 0) (r) = r e cxp (-\utt 2 - hr^j (42) 

Deformed Coulomb potential 

Our last example is the deformed Coulomb potential which can be obtained by introducing 

In this case the equations take the forms 

d 2 (£+!)(£ + 2) 



+ V+- lu t (2£ + 1) + — + e 2 f + (r) = (44a) 



dr 2 r 2 + 4(<? +1)2 



d 2 £(£+!) 



+ V-- u T {2£ + 3) + — — + e 2 /_(r) = (44b) 



<9r 2 r 2 v ' 4(^+1) 
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F± = -7 + ^ 2 + ^ + I(iTIF (45) 

Its groundstate eigenfunction are given by 

The non-relativistic two dimensional Hamiltonian described Coulomb interaction between charged particles, i. e. 
the interaction between conduction electron and donor impurity center when a constant magnetic field is applied 
perpendicular to the plane of motion has been discussed in the literature [H El HI HI . We have obtained relativistic 
interaction and it can be solved quasi-exactly. Recently its spectrum and wavefunction has been computed numerically 



CONCLUSION 



We have obtained analytical solutions of the (2 + 1)- dimensional Dirac equation for a set of potentials in two 
dimensions with the hope that they could be useful in the low dimensional field theory and condensed matter physics. 
The potentials for Dirac equation have been obtained by extending the notion of the Dirac oscillator. In a similar 
manner one can construct the exactly and QES Dirac equation including hyperbolic and trigonometric potentials. It 
has been shown that the (2 + l)-dimcnsional Dirac equation can be transformed in the form of a Schrodinger-likc 
equation and for all exactly and QES Schrodinger equations one can find potentials of the Dirac equation which are 
also exactly solvable or QES. 

Therefore, besides its importance as a new treatment of the construction of the various potentials for Dirac equation, 
the potentials obtained here might be relevant to model some physical problems. Before ending this work a remark 
is in order. It is expected that the models presented here to construct Dirac equation including various potentials 
may provide a good starting point for the study of more realistic models for the low dimensional of field theory and 
condensed matter physics. 
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